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On the origin of W UMa type contact binaries - a new method for 
computation of initial masses 
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ABSTRACT 

W UMa type binaries have two defining characteristics. These are (i) the effective temper- 
atures of both components are very similar, and (ii) the secondary (currently less massive) 
component is overluminous for its current mass. We consider the latter to be an indication 
of its mass before the mass-transfer event. For these stars, we define a mass difference (SM) 
between the mass determined from its luminosity and the present mass determined from fit- 
ting the binary orbit. We compare the observed values of the mass difference to stellar models 
with mass-loss. The range of initial secondary masses that we find for observed W UMa type 
binaries is 1.3 — 2.6M . We discover that the A- and the W-subtype contact binaries have 
different ranges of initial secondary masses. Binary systems with an initial mass higher than 
1.8 ± O.IMq become A-subtype while systems with initial masses lower than this become 
W-subtype. Only 6 per cent of systems violate this behavior. We also obtain the initial masses 
of the primaries using the following constraint for the reciprocal of the initial mass ratio: 
< \/q% < 1. The range of initial masses we find for the primaries is 0.2 — 1.5M©, except 
for two systems. Finally in comparing our models to observed systems, we find evidence that 
the mass transfer process is not conservative. We find that only 34 per cent of the mass from 
the secondary is transferred to the primary. The remainder is lost from the system. 

Key words: binaries: close-binaries: eclipsing-stars: evolution-stars: interior-stars: late-type. 
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1 INTRODUCTION 

The evolution and structure of a star are primarily determined by 
its initial mass. Therefore, the structure of stars in contact binaries 
(CBs) can be very different from single stars because of mass trans- 
fer. The component stars of CBs evolve as single stars with initial 
mass (Mi) until the binary interaction occurs. When the mass trans- 
fer starts, the initially more massive component (M21, the current 
secondary) loses mass that is transferred to the initially less mas- 
sive component (Mii, the current primary) and the rest leaves the 
system. Because of the interaction, the mass ratio of the system 
reverses. After mass transfer stops or sufficiently slows down, the 
secondary stars are found to be overluminous for a main-sequence 
(MS) star of the same mass. The reason for this luminosity excess 
is that the luminosity of such stars depends on their initial mass as 
much as their present mass. In this study, we develop a method to 
find initial masses of the secondary stars of W UMa type binaries 
in terms of the luminosity excess. 

W UMa type CBs are among the best studied stars in the liter- 
ature. For about 100 CBs (51 A-subtypes and 49 W-subtypes), ac- 
curate dimensions of component stars from photometric and spec- 
troscopic data are successfully derived by different investigators. 
The subtypes are defined by their light-curves (Binnendijk 1970). 
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We compute initial masses of secondary and primary components 
of these systems. Our sample include s well-determined parame - 
ters of CBs from the co mpil ations of Yaku t & Eggletor 1 J2005b . 



Gazeas et alj < |2005|.|2006|) and lZola et alj d2005l 12010) 

Beside the luminosity excess of the secondary components 
in the W UMa type binaries, another very important feature of 
these systems is that the effective temperatures, T c , of compo- 
nents are almost identical. However, the slight difference between 
effective temperatures of components is systematic between the 
two subtypes. For the A-subtype systems, the effective temperature 
of the massive component (Ti e ) is higher than the effective tem- 
perature of the lower mass component (Xbe); the opposite is true 
for the W-subtype systems. There is debate in the literature con- 
cerning whether the subty p es represent an evolutiona r y sequence 
dAwadalla & Hannd l2005k iGazeas & Niarchosl IgOOfj : lEker et alj 
|2006|) . According to lMaceroni & van't Veerl l 19961) . however, there 
is no evolutionary link between most of the A-subtype systems and 
the W-subtype CBs. In this paper, we aim to confirm that the sub- 
types are determined by their initial parameters. If initial parame- 
ters of A- and W-subtypes are completely different, then they are 
not an evolutionary sequence. 

There are two evolutionary pathways that lead to a star filling 
its Roche lobe. These are nuclear an d angular momentum evolution 
dHilditch. King & McFarianelll989h . It is thought that angular mo- 
mentu m evolution via magnetic braking dSchatzmanll 19621 : iMestell 
Il968l) is the primary mechanism for the formation and evolution of 
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WUMa type CBs dOkamoto & Satalll97Cl : lvan't VeeJI 197g| : IVilhul 
Il982h . However, both these two mechanisms affect the kind of W 
UMa binaries that are formed. Also because nuclear and angular 
momentum evolution are very strongly dependent on stellar mass, 
one may expect that the initial masses of components of W UMa 
type binaries are the key parameters for our understanding of the 
formation and evolution of W UMa type binaries. 

Each component of a CB fills its Roche lobe and is under 
the gravitational interaction of the other. Apart from this, the more 
massive primary component is an MS star and the luminosity of 
the secondary component is higher than the expected MS luminos- 
ity dYang & Lij|200ltlWebbinMl2003l : lLTet al.ll2008l) . This excess 
may be due to large amount of energy tran sfer that occu rs (see Sec- 
tion 2) from the prim ary to the secondary dLucvlll968allbT IWebbinld 
Li et alj|2008l) . Ho wever, the evolutionary status of CB s is 
still unclear. According to lHilditch. King & McFarlanel dl988h . the 
primary components of the W-subtype CBs are not evolved MS 
stars, and the A-subtype CBs are close to the terminal-age of the 
main sequence (TAMS). 

Although there are many studies in the literature about con- 
tact systems, the contact models are still not able to satisfy all 
the observational constraints. It is still an open question whether 
CBs are in thermal equilibrium or not. It is generally accepted 
that a contact configuration in thermal equilibrium is not possi- 
ble and therefore CBs suffer from thermal r elaxation oscillation 
jFlannervlll974 lOianl l200ll: lKahlerll2002allbl |2004 iKreiner et alj 
20031 : IWebbinMl2003l: lli. Han & Zhand 12004 : 1 Yakut & Eggleton] 
20051 : IPaczvnski et al.ll200fj) . However, according to some authors, 
a contact configuration can be in thermal equilibrium. They argue 
that contact systems can achieve thermal equilibrium if the sec- 
ondary, currently less m assive compone n t, is more ev olved than 
the massive component dStepieftoOoiH, 120091 l201lh . In such a 
model, CB stars originate from detached close binaries with initial 
orbital periods close to a few days and the system was exposed to 
mass exchange with a mass ratio reversal in the past. According to 
these assumptions, energy exchange does not influ e nce stellar rad ii 
and occurs in thermal equilibrium fKahler 2004; Stepieri 1201 ll ). 
Large-scale energy transfer from the primary to the se condary com - 
ponents solves several outstanding problems in CBs dLucvlll968ah . 
However, the problem of energy transfer requires an intrinsic driv- 
ing mechanism and its effect on the evolution of the component 
stars and the system is unknown. These processes are very poorly 
understood because of the complexity of detailed hydrodynamic 
treatment of the flow and also due to the unknown driving mecha- 
nism. 

This paper is organized as follows. We outline the physical 
properties of CBs in Section 2 and we discuss the basic properties 
of CBs in the mass-luminosity (M — L) and the mass-radius (M — 
R) diagrams. In Section 3, we develop and apply a comprehensive 
method for computation of initial masses of CBs based on stellar 
modelling with mass-loss. In Section 4, we apply this method to 
the CBs and discuss the results. Finally, in Section 5, we draw our 
conclusions. 



2 THE PHYSICAL PROPERTIES OF CONTACT 
BINARIES AND THE PROBLEM OF ENERGY 
TRANSFER 

The most reliable accurate parameters (mass, radius and luminos- 
ity) of binary stars are computed by using both photometric and 
spectroscopic data. We compiled 100 such well known CBs from 



the literature. Their basic properties are listed in Table A 1 of the 
Appendix. These properties give clues to evolutionary phases of 
the primary components, at least. The primary components of CBs 
in M — L and M — R diagrams are much closer to that of the 
normal MS stars than the secondary components. Their luminosi- 
ties and radii are plotted with respect to their masses in Figs 1(a) 
and 1(b), respectively. For comparison, the components of the well- 
determined de tached eclipsing binaries (+; herea fter DEBs) are 
also included dTorres. Andersen & Gimenez 2010:). The solid line 
in these figures shows zero-age M S (ZAMS), derived from models 
constructed with the MESA code dPaxton et al .11201 lT) . There is an 
agreement between the luminosities of the primaries of CBs and 
that of the components of DEBs. If we compare the radii of the pri- 
mary stars of CBs in Fig. 1(b), we can confirm that these two groups 
of stars are in very good agreement with the DEB stars. Fig. 1(b) 
also shows the TAMS line (as a dashed line). Most of the primaries 
of CBs occur between the ZAMS and the TAMS lines. 

The M — L and M — R diagrams of the secondary components 
of W UMa type CBs are plotted in Figs 1(c) and 1(d), respectively. 
Unlike the primaries shown in Figs 1(a) and Figs 1(b), there is no 
agreement between the secondary components and the components 
of DEBs in M — L and M — R diagrams. Therefore, the secondary 
components of W UMa type CBs must have a very different struc- 
ture or evolutionary path. All the secondary components are above 
the ZAMS line in both Figs 1(c) and 1(d). The mean masses of 
the secondary components of the A- and W-subtype CBs are 0.38 
and 0.47Mq , respectively. Despite these low masses, the secondary 
stars are very bright and very large in comparison to their MS coun- 
terparts of the same mass. This excess in luminosity and radius may 
be explicable by the early evolution of the systems. 

In the literature, the close T c s of components of CBs have 
been explained as a result of an e nergy transfer from the primary 
to the secondary star dSmitbll984l) . According to this approach, the 
luminosity excess of the secondary components (8L2 = Li — M| ) 
must compensate the lack of luminosity of the primary components 
(<5Li = Li- Mi). If this is the case, there must be a correlation 
between SLi and 8L2 ■ In Fig. 2, 8L2 is plotted with respect to 8L\ . 
It should be noted that 8L2 is always greater than zero and SLi is 
in general (four fifths of systems) less than zero. If the luminosity 
excess of the secondary stars compensates the lack of luminosity of 
the primary stars, the position of the CB components should be on 
the 8L2 = — SLi line (solid line). There is a large scatter for the 
A-subtype CBs (circles) and therefore no evidence for the energy 
transfer process is shown for this subtype. For the W-subtype CBs 
(filled circles), however, some of the systems with small values of 
8L1 and 8L2 form a sequence on the line 8L2 = —8L\. The se- 
quence is very clear from the origin to the point with 8L2 = ILq 
and SLi — —ILq. The details of this range are also shown in Fig. 
2. 

The energy transfer process does not work, at least, for the 
A-subtype CBs. For the W-subtype CBs, however, there are some 
systems (about 12) in which the transfer process may play a part 
for very close effective temperatures. For some binaries, however, 
both components have excess luminosity which the energy transfer 
process cannot explain. 



3 METHOD FOR COMPUTATION OF INITIAL MASS 
FROM LUMINOSITY EXCESS 

As one of the essential parameters of a star, luminosity depends 
on the nuclear burning in the core. For a single star, this is deter- 
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Figure 1. Mass-luminosity and mass-radius relation for the primaries of the A- (circle) and the W- subtype (filled circle) CBs are plotted in panels (a) and (b), 
respectively. For comparison, the components of the well-known DEBs (+) are also plotted. Solar units are used. A solid line indicates the ZAMS line and 
a dashed line the TAMS line, obtained by using the MESA stellar evolution code. Panels (c) and (d) are same as (a) and (b), respectively, but for secondary 
components of W UMa type CBs. 



mined by the initial mass of the star. For a component star which 
has experienced mass transfer (and loss), 

M = Mi + AM, (1) 

where AAf is the mass lost or transferred from the secondary, i.e. 
the mass difference between the initial (Mi) and the final (Mf) 
masses. The physical conditions in the central regions of the sec- 
ondary star, which yield luminosity higher than Mf, are similar to 
that of a star with a mass lower than Mi but higher than Mf. More 
precisely, the central physical conditions of the secondary star are 
very close to a single star of mass Ml ~ L ' 25 (in solar units). For 
Ml, which is the mass of an isolated star with the same luminosity, 
one can use a more realistic mass-luminosity relation. The M — L 
relation we employ is 

L — 1. 49 M 4 216 , (2) 

derived from the grids obtained bv lYildizI d2012l) for TAMS mod- 
els. The value of mass according to luminosity strongly depends 
on the evolutionary phase of a mass losing star at the beginning 
of the mass loss/transfer process. For the computations relating to 



the secondary stars of CBs, for now, we assume that mass transfer 
starts near the TAMS phase of the old massive component. There- 
fore, equation (2) is used throughout this paper [see Section 5.3 for 
correction due to deviation from equation (2)]. For a typical TAMS 
star, the quantity M L defined by Ml = (£ 2 /1.49) 1/4 - 216 for the 
secondary star computed from equation (2) is equal to the observed 
M2 = Mf. If there is a luminosity excess, then Ml > M%. We 
define the mass difference between these two masses as 5M: 

5M = M L -M 2 . (3) 

According to our assumption, the luminosity excess of the sec- 
ondary component of a CB is related to its initial mass as much as 
to its present mass, and therefore to the amount of total mass lost by 
the secondary (AM). This means that there must be a correlation 
between AM and SM: AM = f(SM). Then the initial mass of 
the secondary component of a CB is given as 

M 2i = Mi + f(SM) (4) 

We note that M 2 i is a function only of the observable quantities L2 
and M%. The masses, luminosities and radii are in solar units. 
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Figure 2. Excess in luminosity of the secondary component (8L2 = L2 — 
Mi) is plotted with respect to luminosity lack (or excess) of the primary 
components (5L± = L\ — Mf). If the energy transfer is the dominant 
operating mechanism, all the CBs should make a sequence on the SL2 = 
—SLi line (solid line). This is the case for only some of W-subtype CBs. 



If one finds Mm, then Mij can be computed by using the con- 
straint on the reciprocal of initial mass ratio 1 /gi : < 1/g; < 1. 
Mii is a function of AAI — M21 — M2 and mass lost by the system 
(Afi os t). The increase in the primary mass can be written as 

Mi - Mu = AM - M ost = AM(1 - 7). (5) 

where 7 is the ratio Mj ost /AM. 7 = corresponds to the con- 
servative case and 7 = 1 is the case in which all the mass lost by 
secondary component is lost by the system. The reciprocal of the 
initial mass ratio as a function of current masses, M21, and 7 is 
given as 

(M 2i -M 2 )(l- 7 ) 



t , Mi; Mi 
1/9i= M^ = — 



M 2i 



(6) 



The unknown value of 7 can be found from the constraint < 
1/gi < 1 (see section 4.2). Our next task is to find an expression 
for/(<5M). 

3.1 Method derived from the mesa models 

3.1.1 Initial masses of the secondary components 

As discussed above, the secondary components of CBs have a com- 
pletely different interior from the normal MS stars because their 
initial mass is significantly higher than the present mass. This sit- 
uation forces the secondary components to have a different evolu- 
tionary path from single MS stars without mass-loss. Change in 
mass essentially causes a breakdown in hydrostatic equilibrium, 
which is then restored by expansion of the central regions. We con- 
struct models in order to assess the ultimate effect of change in 
mass on the structure of stars by using the MESA stellar evolution 
code. From these models, we aim to find an expression for AM as 
a function of S AI. 

Models with constant and variable masses are constructed by 
using the MESA co de. The details of this code can be found in 
IPaxton et al.l fcOllh . The chemical composition and convective pa- 
rameter for all the models are taken as X=0.700, Z=0.02, a=1.89. 
Stellar models are constr ucted with the solar mixture given by 
iGrevesse & Sauval dl998l) . The radiative opacity is derived from 



Table 1. Model results with mass-loss to obtain an expression for the 
difference between initial and final masses. The final mass is chosen as 
0.5Mq and the mass transfer begins when central hydrogen abundance 

x c = icr 6 . 



Model 


Mi 


M f 




M L 


SM 


AM 


No 


(M ) 


(M©) 


(L©) 


(Mq) 


(M©) 


(M ) 


1 


1.20 


0.50 


0.251 


0.708 


0.208 


0.700 


2 


1.50 


0.50 


0.359 


0.774 


0.274 


1.000 


3 


1.80 


0.50 


0.638 


0.894 


0.394 


1.300 


4 


2.50 


0.50 


2.884 


1.303 


0.803 


2.000 


5 


3.00 


0.50 


6.457 


1.594 


1.094 


2.500 


6 


3.50 


0.50 


11.220 


1.830 


1.330 


3.000 


. Same as Table 1 , but for the final 


mass Mf = 0.8M Q . 


Model 


Mi 


Mf 




M L 


SM 


AM 


No 


(M ) 


(M Q ) 


(L ) 


(M Q ) 


(M ) 


(Mq) 


1 


1.20 


0.80 


0.893 


0.972 


0.172 


0.400 


2 


1.50 


0.80 


1.064 


1.016 


0.216 


0.700 


3 


1.80 


0.80 


1.346 


1.077 


0.277 


1.000 


4 


2.50 


0.80 


3.926 


1.408 


0.608 


1.700 


5 


3.00 


0.80 


8.630 


1.714 


0.914 


2.200 


6 


3.50 


0.80 


17.378 


2.042 


1.242 


2.700 



recent OPAL tables dlglesias & Rogers1ll996t). im plemented by the 
low-temperature tables of Ferguson et al.l (120051). The the r monu- 
clea r reactions rates are compu ted by using lAngulo et al.l dl999h 
and ICaughlan & Fowleil dl988l) . a nd the standard m ixing length 
theory is employed for convection dCox & G iuli 1968). The MESA 
p — T tables are based on th e 2005 update of the OPAL EOS tables 
dRogers & Navfonovll2002h . 

The progenitors of the secondary stars of CBs should be 
evolved stars (Hilditch et al. 1988), otherwise an excess in lumi- 
nosity cannot be observed. Therefore, when we construct our stel- 
lar models we assume that mass transfer begins near the TAMS 
which we take to be when the central hydrogen abundance is 
10 -6 . The basic properties of the constructed models with differ- 
ent initial masses are listed in Table 1. The range of initial masses 
is 1.2 — 3.5Mq. For all of these models, the mass loss rate is 
lO _8 M0yr~ 1 and the final mass is 0.5M©. After several tries for 
different mass-loss rates, we deduce that the results are independent 
of the mass-loss rate. The corresponding mass to the luminosity 
(Ml) after mass transfer is computed from equation (2). The mass 
difference SA1 is as given in equation (3) and AAI — Mi — Mi. 
We also construct models with Mf = O.8M0. The basic properties 
of these models are given in Table 2. 

The initial mass is plotted with respect to Ml in Fig. 3 for 
Me = O.5M and 0.8M©. The relation between the initial and lu- 
minosity masses is almost linear. In both sets of models with either 
M f = 0.5M© or O.8M , we note that the higher the AM, the 
higher is the 5AI. In Fig. 4, A/A/ is plotted with respect to 8M (in 
units of M Q ). Both sets of models with Mf = O.5M (filled circle) 
and Mf = O.8M (circle) have the same relation between AM and 
SM. The solid line shows the fitted line [2.50(5M - 0.07) o e4 ] for 
both data sets. It is remarkable to confirm that the relation is not 
a linear one. Also shown in Fig. 4 is the dotted line for the linear 
relationship AM = 2.05M. 

The initial mass of the secondary stars of CBs is then com- 
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Figure 3. The relation between Ml and Mi for the final masses 0.5 (filled 
circle) and 0.8 (circle) Mq. 
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Figure 4. AM as a function of 8M for the MESA models with final masses 
M t = 0.5 (filled circle) and Mf = O.8M (circle). The solid line repre- 
sents the fitting curve AM = 2.50(<5M - 0.07) 64 and dotted line is for 
AM = 2.0SM. 



puted by using 

M 2i = M 2 + AM = M 2 + 2.50(JM - 0.07) ' 64 (7) 
where the masses are in units of solar mass. 



3.1.2 Uncertainties in the initial masses of the secondary 
components 

The initial masses of the secondary components are computed from 
L 2 and M 2 . Therefore, the uncertainty in M 2 \ (AM 2 i) arises from 
uncertainties AM 2 and AL 2 . Define / = AM = 2.50(5M - 
0.07) ' 64 . Then, / is a function of Ml and M 2 . In terms of these 
quantities and their uncertainties, the uncertainty in M 2 \ can be 
written as 



AM 2 i = AM 2 



df dML -AL 2+ ^ AM2 . 



DM L dL 2 



dM 2 



(8) 



Figure 5. Initial mass of the secondary components of CBs. Circles are 
for A-subtype and the filled circles are for W-subtype. The solid black line 
represents the derived transition mass with the dashed lines showing the lc 
error bars. Apart from a few systems, the A-subtype CBs have M 2 \ higher 
than 1.80 M and the W-subtype CBs have M 2i lower than 1.8OM . 



Using the definition of 8 M given in equation (3), we find that 



df 
8M L 



df 



8M 2 



1.59(<5M - 0.07) ~ 



(9) 



From the employed M — L relation, we obtain 
dM L 0.22 



dL 2 



(10) 



The explicit form of AM 2 i as a function of AM 2 and AL 2 is 
1.59 /0.22AL2 



AM 2i = AM 2 + 



{5M - 0.07) 36 



L0.76 



+ AM; 



(11) 



4 RESULTS AND DISCUSSION 

4.1 The initial masses of the secondary components 

The initial mass of the secondary component of a CB is computed 
using the method from the previous section. The results are listed 
in Table Al. If we compare an A-subtype with a W-subtype CB of 
the same M 2 (present secondary mass), then we find that the initial 
mass of the secondary of an A-subtype is 0.5Mq higher than that 
of W-subtype. This means that the (present) secondary components 
of the A-subtype have lost/transferred much more mass than those 
of the W-subtype. 

In Fig. 5, the initial mass of the secondaries, derived by using 
the method based on MESA models, is plotted with respect to the 
present mass of the secondary components. The initial masses are 
plotted with their uncertainties computed by a method presented in 
the previous section. It can clearly be seen that subtypes are de- 
termined by the initial mass of the secondaries. If M 2 \ > 1.8 ± 
0.1M Q , the system becomes A-subtype. If M 2i < 1.8 ± 0.1M©, 
then the system exhibits the W-subtype properties. There are just 
a few systems that do not obey this prescription (see below). The 
mean initial mass of the secondaries is M 2 \a = 1.97M© for the 
A- and JI 2iw = 1.56M© for the W-subtype CBs. The underly- 
ing mechanism for this separation must be based on competition 
between the orbital and the nuclear time-scales. 
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There are only two W-subtype systems (Wa) having an initial 
mass higher than M 2 iA- They are V728 Her and V402 Aur. These 
binaries are the hottest of the W-subtype CBs. Their M — L re- 
lation and angular momentum properties are similar to that of the 
A-subtype rather than the W-subtype. The number of the A-subtype 
binaries (Aw) having an initial mass significantly less than M 2 iW, 
however, is four. These binaries are EQ Tau, 00 Aql, V508 Oph 
and TZ Boo. They have T c about 5600 K [log(T c ) = 3.75] and 
are the coolest of the A-subtype CBs. The O'Connell effect may be 
responsible for this atypical situation if it is not due to uncertain- 
ties in the observed T c s. However, in addition to the initial mass of 
the secondary components, another secondary effect may also play 
a role in the formation of W UMa subtypes. For example, initial 
orbital angular momentum may be important. 

Putting the atypical CBs aside, the initial mass range for the 
secondary components is 1.7 — 2.6M© for A-subtype and 0.7 — 
1.9M© for W-subtype. However, for most of the CBs the range is 
1.3 - 2.2M . 

CC Com is the system with the minimum secondary initial 
mass, 0.86 Mq. The nuclear time-scale of such a star is about 
16 Gyr, greater than the age of the galaxy. Therefore, this system 
(or most of the W-subtype CBs) must have some specific feature, 
which makes its fundamental properties unusual. As a matter of 
fact, both components of CC Com are brighter than ZAMS stars 
and fainter than TAMS stars. Formation of binaries with M 2 i < 
1.25Mq seems to be different from the other CBs. They have ei- 
ther a different value of 7 or different structure for the progenitors 
of the secondary components during the mass transfer process. For 
example, the progenitors may not be as evolved (not TAMS) as we 
assume (see Section 5.3). 

The most realistic uncertainties in M 2 i are computed from 
equation (11) and given in Table Al. For some of the binaries, for 
example, HV UMa, the uncertainty AM 2 i is very high due to high 
uncertainties in M 2 and L 2 . However, this result does not change 
the situation much because the uncertainties are within an accept- 
able range for most of the binaries in our sample. 

4.2 Initial masses of the primary components 

The reciprocal mass ratio 1 /qi is the physical constraint for the de- 
termination of initial masses of the primary components. The fitting 
parameter 7 (equation 5) is the mass lost by the system in units of 
AM. First, we consider A-subtype CBs. If 7 has a small value, 
then some A-subtype systems have an unreasonably negative value 
for Mn. In order to keep Mu > (1/qi > 0) for all A-subtype 
CBs, the minimum value of 7 is 0.50. The maximum value of 7 
for 1/qi < 1, however, is 0.87. These two values of 7 give mean 
value for 1/qi as 0.368 and 0.666, respectively. For 1 /g i = 0.5, 
7 = 0.664. This value of 7 yields Mu ranging from 0.3 to 1.5M© 
for A-subtypes. This means that the mass lost by these subtypes 
is higher than the transferred mass: 66 per cent of mass lost by 
the secondary components is lost by the system and 34 per cent is 
transferred to the primary components. 

1/qi of A-subtype CBs with 7 = 0.664 is plotted with re- 
spect to M2 in Fig. 6. The reciprocal of the initial mass ratio is 
well within the expected range. For 1 /qi of W-subtype CBs, we use 
the same value of 7. On average, the A-subtype CBs have higher 
1/qi than the W-subtype CBs. There are only 3 violating systems. 
These binaries are the W-subtype CBs (see below) and have very 
low 8M values, less than 0.2M©. They may not be TAMS stars but 
are somehow in between the ZAMS and the TAMS lines. 

The value of 7 we found is a mean value. Indeed, there is no 
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Figure 6. l/qi computed from the expressions derived from the models by 
using MESA code is plotted with respect to M 2 . All the A-subtype systems 
are in the range < 1/qi < 1; just three W-subtype systems are greater 
than one. 



restriction for a single value of 7 for all CBs. However, what is 
important is that there is a solution with a single value of 7 for 
almost all CBs, and this is certainly true for the A-subtype CBs. It 
may be concluded that the value of 7 does not dramatically change 
much from system to system. 

The initial masses of the primary components are listed 
in Table Al. The range of primary initial mass of the CBs 
is 0.2 — 1.5Mq, apart from two binaries (HV UMa and 
V376 And). This is the mass range in which magnetized stel- 
lar wind is an efficient mechanism for angular momentum loss 
dTutukov. Dremova & Svechnikovll2004l) . The mean value of the 
initial primary masses is Mn — 0.99Mq for the A- and Mu = 
O.84M for the W-subtype. The uncertainties in Mu are also given 
in Table Al. AMu is taken as AM 2 i/qi. 



4.3 Correction of the initial masses of components of CBs 
with SM < 0.35M o 

As previously stated, the initial mass of some secondaries is very 
small for some CBs, for example CC Com. In Fig. 7, the recip- 
rocal of initial mass ratio 1/qi is plotted with respect to 8 M both 
for the A- and the W-subtype CBs. For A-subtype, 1/gi does not 
show any dependence on 5M. For W-subtype, however, 1/qi is de- 
pendent on 5M if 5M < 0.4Mq. The solid line shows the curve 
fit (-3.055M - 25 + 2.80) for this range. This relation reflects our 
assumption about the evolutionary phase of secondary components 
at the beginning of mass transfer. It seems that for some of the 
W-subtype binaries the secondary components are not close to the 
TAMS line. Therefore, equation (7) is not valid for these CBs. Us- 
ing the curve (-3.055M 25 + 2.80) of Fig. 7, we obtain a new 
equation for M 2l , and Mu is computed as in Section 5.2. This 
equation removes the slope in Fig. 7. 

The corrected initial masses of the secondary and the primary 
components of these CBs are listed in Table 3. The upper limit for 
the initial mass of the secondaries is 1.90Mq, consistent with the 
results given above. The minimum initial mass of the secondaries 
is with small SM pertaining to BH Cas, 1.36Mq. 
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Figure 7. 1/qi computed from the expressions derived from the models by 
using the MESA code is plotted with respect to SM. l/q\ for the A-subtype 
CBs has a very uniform distribution. For the W-subtype CBs, however, this 
is not the case. For small values of SM, 1/qi is dependent on SM, The solid 
line shows the fitted curve for the W-subtype CBs (-3.055M - 25 + 2.80). 



Table 3. The corrected masses for the secondary and the primary masses 
of the W-subtype CBs with SM < 0.35Mq. The correction is com- 
puted according to the expression obtained from the curve in Fig. 7 
(-3.055M 025 + 2.80). 



Star 


Mi 


M 2 


M 2i 


M u 


l/9i 


*ms (Gyr) 


DN Cam 


1.85 


0.82 


1.86 


1.50 


0.81 


1.25 


EF Boo 


1.61 


0.82 


1.90 


1.25 


0.66 


1.17 


ET Leo 


1.59 


0.54 


1.55 


1.25 


0.81 


2.19 


AAUMa 


1.56 


0.85 


1.88 


1.21 


0.65 


1.20 


YY Eri 


1.55 


0.62 


1.53 


1.24 


0.81 


2.28 


EROri 


1.53 


0.98 


1.74 


1.27 


0.73 


1.52 


BB Peg 


1.42 


0.55 


1.59 


1.07 


0.67 


2.03 


AEPhe 


1.38 


0.63 


1.69 


1.02 


0.61 


1.68 


WUMa 


1.35 


0.69 


1.77 


0.99 


0.56 


1.44 


AM Leo 


1.29 


0.59 


1.64 


0.94 


0.57 


1.85 


V781 Tau 


1.29 


0.57 


1.61 


0.94 


0.59 


1.96 


RZCom 


1.23 


0.55 


1.58 


0.89 


0.56 


2.09 


AO Cam 


1.12 


0.49 


1.54 


0.77 


0.50 


2.26 


GZAnd 


1.12 


0.59 


1.68 


0.75 


0.45 


1.69 


TWCet 


1.06 


0.61 


1.74 


0.68 


0.39 


1.52 


BX Peg 


1.02 


0.38 


1.42 


0.67 


0.47 


2.90 


AB And 


1.01 


0.49 


1.54 


0.66 


0.43 


2.25 


SWLac 


0.98 


0.78 


1.88 


0.61 


0.33 


1.21 


VWCep 


0.93 


0.40 


1.48 


0.57 


0.38 


2.56 


V757 Cen 


0.88 


0.59 


1.74 


0.49 


0.28 


1.52 


CC Com 


0.72 


0.38 


1.49 


0.35 


0.23 


2.50 


XYLeo 


0.76 


0.46 


1.68 


0.35 


0.21 


1.69 


V523 Cas 


0.75 


0.38 


1.53 


0.36 


0.24 


2.28 


BH Cas 


0.74 


0.35 


1.36 


0.40 


0.29 


3.32 


RW Dor 


0.64 


0.43 


1.63 


0.24 


0.14 


1.86 



4.4 Some notes on evolutionary time-scales and cluster 
member CBs 

As mentioned above, the ultimate fate of a binary depends on its 
total mass and how this mass is shared between its components. If 
mass of the massive component (M 2 i) is high enough, then nuclear 
evolution is very fast compared to orbital evolution. Therefore, the 
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Table 4. Cluster member CBs. 



Star 




M 2i 


tMs(Gyr) 


id (Gyr) 


P(day) 


Subtype 


Cluster 


TXCnc 


0.50 


1.70 


1.6 


0.7 


0.3830 


W 


Praesepe 


AHVir 


0.94 


1.67 


1.7 


3.0 


0.4075 


W 


Wolf 630 


QXAnd 


0.62 


1.92 


1.1 


2.0 


0.4118 


A 


NGC 752 


AHCnc 


0.93 


1.85 


1.3 


4.0 


0.3604 


A 


M67 



Table S. Mass intervals for final products. 



minM 2 i(M Q ) 


maxM 2i (M ) 


Mu (M ) 


tmm(Gyr) 


Product 


1.3±0.1 


1.8±0.1 


0.2-1.5 


1.4 


W-subtype 


1.8±0.1 


2.6±0.1 


0.2-1.5 


0.5-1.4 


A-subtype 



lifetime of this star is so short that the system cannot be a CB . In 
this respect, a very interesting result we obtain is that there is an up- 
per limit for the mass of the massive components of CBs, 2.61Mq 
(V921 Her). This result leads us to consider binaries with compo- 
nents having a mass higher than 2.6IM0 as candidate systems for 
binaries with compact objects, for example. 

Ages of the W UMa type CBs ar e investigated by many dif- 
ferent groups us i ng different methods ( Guinan & Brads tree t il 19881: 
van't Veerlll996l : iBilir at al.ll20ol ILL Han & Zhangll2005l: Li et al.l 
20071) and range from 0.1 - 1.0 to 7.2Gyr. The life-time of a W 
UMa type CB can roughly be divided into three different phases. 
These are detached, semi-detached and contact phases. From the 
initial masses, at least for the A-subtype CBs, we find how long 
the detached phase is: tdet — tins- The MS age for a star with 
2.61Mq is 0.46Gyr. This time is likely to be very short for angular 
momentum evolution without mass-transfer. The MS age of a star 
with transition mass (1.8Mq ) from the A- to the W-subtype CBs is 
1.37Gyr (EX Leo). For such systems, the nuclear and angular mo- 
mentum time-scales are comparable. For some W-subtype systems, 
the nuclear time scale is liable to be so long that the components 
evolve into contact phase before significant nuclear evolution of the 
massive component. The MS lifetime for BH Cas, for example is 
3.32Gyr. Such a time is very long for angular momentum evolution 
of BH Cas and similar CBs with two late-type components. 

According to the results, we deduce that ages of the A-subtype 
binaries (t„) are greater than the MS lifetime of the progenitor of 
their secondaries (4ms )• F° r the W-subtype binaries, however, t s 
may also be less than £ms because it seems that the mass transfer 
for some CBs starts before the secondary component (M 2 i) reaches 
the TAMS line. In order to test this point, we consider CBs in clus- 
ters. There are four such binaries in our CB sample. Two of them 
are the W-subtype CBs. The W-subtype binaries TX Cnc and AH 
Vir are members of Praesepe and Wolf 630, respectively. The A- 
subtype CBs are QX And of NGC 752 and AH Cnc of M67. The 
MS lifetime according to M 2 i (£ms) of the binaries and age of 
clusters are listed in Table 4, respectively. While £ms is less than 
ts = icius (cluster age) for A-subtype, tMS for the W-subtype bi- 
naries is less than t B for AH Vir and greater than for TX Cnc. Al- 
though there are only four well-known cluster member systems, 
these results are consistent with our expectations. 

Table 5 summarizes the results obtained for the initial masses 
of the A- and the W-subtype CBs. The MS lifetime according to 
M 2 i is given in the fourth column. The age of the A-subtypes must 
be greater than this age. For the W-subtype CBs, however, age can 
be either greater or smaller than the given value 1.4Gyr, depending 
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5 CONCLUSION 

In the present study we have investigated the physical properties of 
CBs. The approximate initial masses of the secondary and primary 
components are computed, for the first time ever, by using the basic 
parameters of 100 W UMa type CBs (51 A- and 49 W-subtype). 
The initial mass of the secondary is computed from its luminosity 
excess by a method based on the expression derived from stellar 
models with mass-loss. We have also deduced that only one-third 
of the mass lost by the secondary components is transferred to the 
primary components. The remaining part leaves the systems. The 
range of initial masses of the secondaries of CBs is 1.3 — 2.6M0. 
The upper part of this range is occupied by the A-subtype CBs and 
the W-subtype CBs are in the lower part. The transition occurs at 
about 1.8±O.1M . 

The main assumption in our method is that the mass trans- 
fer starts near or after the TAMS phase of the massive component 
(progenitor of the secondary component). This assumption seems 
fairly good for the A-subtypes. For the W-subtype, however, a cor- 
rection is needed. From SAI and the reciprocal mass-ratio diagram 
(Fig. 7), we develop a method and obtain corrected initial masses 
for the case of small SAL The corrected initial masses of the sec- 
ondaries range from 1.3 — 1.9M0. One of the important outcomes 
of the present study is that there is also a lower limit for the initial 
mass of the primary components. The reason for this result should 
be studied in the context of angular momentum evolution versus 
nuclear evolution. 

For the binary systems with an initial mass higher than 
1.8Mq, although the angular momentum is lost by the less mas- 
sive component the nuclear evolution is principally the mechanism 
responsible for the Roche lobe filling process. Therefore, semi- 
detached systems wifhAfe = 1.8 — 2.61Mq must yield A-subtype 
W UMa CBs. The binary systems with M21 less than 1.8Mq evolve 
into the contact phase due to the rapid angular momentum evolu- 
tion. The Roche lobes of their components contract so that they 
form the W-subtype contact or near-contact binaries despite the 
components not being much evolved in the MS phase. 

The primary components of W UMa type CBs have an initial 
mass range from 0.2 to 1.5Mq . A star with mass in this range loses 
its spin angular momentum relatively rapidly. If it has a tidally in- 
teracting companion, orbital angular momentum is lost by the sys- 
tem. The precursors of W UMa type CBs are the systems in which 
nuclear evolution of the massive components is in competition with 
the angular momentum loss from the less massive components. In 
the precursors of the W-subtype W UMa CBs, both components are 
effective in the angular momentum loss process. 

The method we develop for computation of the initial masses 
of the secondary components is highly suitable for a complete er- 
ror analysis. From the derivatives of the function derived from the 
fitting curve for AM as a function of L2 and Ma (equation 8), we 
successfully compute AA'hi in terms of the observed uncertainties 
AL2 and AM2 (equation 11). 

We also consider cluster member W UMa type CBs. Although 
the number of CBs with accurate dimensions available in the liter- 
ature is very small (4 systems), the results support our findings on 
initial masses and time-scales for the A- and W-subtype CBs. 
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APPENDIX A: BASIC PROPERTIES OF W UMA TYPE BINARIES 



Table Al: Fundamental properties of the A- and W-subtype CBs. Columns are organized 
as name, subtype, period, mass ratio, effective temperatures of primary and secondary, 
primary and secondary radii, luminosities of primary and secondary, primary and secondary 
masses, initial masses of primary and secondary before the mass transfer (see Section 5), 
reciprocal of initial mass ratio, MS life time of secondary component according to initial 
mass and numbers for references. 



Star 

HVUMa 
V376 And 
RRCen 
V921 Her 
UZ Leo 
V535 Ara 
AQTuc 
EFDra 
V2388 Oph 
AWUMa 
V776 Cas 
V592 Per 
XZ Leo 
DKCyg 
NNVir 
»?CrA 
RZTau 
V401 Cyg 
AQPsc 
AH Aur 
V839 Oph 
FPBoo 
VI 073 Cyg 
EX Leo 
AHCnc 
APLeo 
FGHya 
UXEri 
YY CrB 
CKBoo 
V566 Oph 
TYC1174 a 
CN Hyi 
GR Vir 
DZ Psc 
RTLMi 
HNUma 
V410 Aur 
SXCrv 
EQTau 
YSex 
V2357 Oph 
V404 Peg 
QXAnd 
VZLib 
00 Aql 
V508 Oph 
DXTuc 
TVMus 
XYBoo 
TZBoo 
DN Cam 
V728 Her 
V402 Aur 
EF Boo 
ET Leo 
AAUMa 



t P 
d 

A0.711 
A 0.799 
A 0.606 
A 0.877 
A 0.618 
A 0.629 
A 0.595 
A 0.424 
A 0.802 
A 0.439 
A 0.440 
A 0.716 
A 0.488 
A 0.471 
A 0.481 
A 0.591 
A 0.416 
A 0.583 
A 0.476 
A 0.494 
A 0.409 
A 0.640 
A 0.786 
A 0.409 
A 0.360 
A 0.430 
A 0.328 
A 0.445 
A 0.377 
A 0.355 
A 0.410 
A 0.389 
A 0.456 
A 0.347 
A 0.366 
A 0.375 
A 0.383 
A 0.366 
A 0.317 
A 0.341 
A 0.420 
A 0.416 
A 0.419 
A 0.412 
A 0.358 
A 0.507 
A 0.345 
A 0.377 
A 0.446 
A 0.371 
A 0.298 
W 0.498 
W 0.471 
W 0.604 
W 0.430 
W 0.347 
W 0.468 



(K) 
70()0± 200 
8350± 120 
6918± 319 
7700± 20 
6980± 15 
8200± 500 

6980 

6000 
6900± 23 

7175 
6700± 90 
6800± 40 
7240± 13 

7500 
6900± 15 

7100 
7300± 8 

6700 
6100± 15 
6215±9 

6650 
6980± 14 
6700± 100 
6340± 14 
6300± 90 
6150± 25 
6200± 20 
6100± 28 
6135± 10 

6200 
7000± 40 
65O0± 36 
6500± 200 
6300± 10 
6210± 12 
6200± 58 
6100± 77 
5890± 22 
6340± 32 
5860± 10 

6210 
5640± 25 

6340 
6500± 150 
5770± 21 
5700± 300 

6000 
6250± 37 
5980± 32 

7200 
5888± 271 
6530± 23 

6622 
6700± 31 

6338 
5112± 28 

5929 



T 2 

(g) 
7000± 200 
7335± 120 
6761± 312 
7346± 20 
6830± 15 
8129± 521 

6860 

6054 
6349± 23 

7022 
6725 ± 90 
6020± 40 
6946± 13 

6700 
6925 ± 15 

6640 
7194± 8 

6650 
6124± 15 
6141±9 

6554 
6456± 14 
6590± 100 
6110± 14 
6275 ± 90 
6250± 25 
6519± 20 
6340± 28 
6142± 10 

6291 
6881±40 
6357± 36 
6360± 260 
6163±9 
6187± 12 
6350± 58 
6082± 77 
5983± 22 
6160± 32 
5851± 10 

6093 
5780± 25 

6154 
6421± 155 
5980± 21 
5680± 300 

5830 
6182± 37 
6088± 32 

7102 
5754± 265 
6700± 23 

6787 
6775 ± 31 

6450 
5500± 28 

5965 



Ri 

(Re) 

2.62± 0.25 
2.66± 0.02 
2.24± 0.10 
2.75± 0.02 
2.29± 0.02 
2.09± 0.03 
2.05± 0.11 

1.70 
2.60± 0.02 
1.87±0.05 
1.82± 0.02 
2.25± 0.03 
1.69± 0.02 

1.71 
1.72± 0.01 
2.20± 0.03 
1.56± 0.07 

1.98 
L75± 0.01 
1.85± 0.02 

1.48 
2.31±0.03 
2.24± 0.02 
1.56± 0.01 
1.40± 0.09 
1.48±0.05 
1.44± 0.01 
1.47± 0.01 
1.43± 0.01 

1.45 
1.47± 0.01 
1.45± 0.02 
1.60± 0.02 
1.42± 0.07 
1.47±0.02 
1.29± 0.02 
1.44± 0.01 
1.44± 0.03 
1.31 ±0.02 
L14± 0.01 
1.50± 0.02 
1.39± 0.02 

1.35 
1.40± 0.01 
1.17± 0.05 
1.38± 0.02 

1.06 
1.20± 0.04 
1.47± 0.17 
1.28± 0.05 
0.96± 0.13 
1.78± 0.02 

1.81 
1.98± 0.02 

1.46 
1.36± 0.01 

1.47 



R2 

(Rq) 
1.18± 0.16 
1.55±0.01 
1.07± 0.05 
1.41± 0.01 
1.39± 0.01 
1.23±0.02 
1.32± 0.07 

0.78 
1.30± 0.01 
0.66± 0.02 
0.75± 0.01 
1.47± 0.02 
1.00± 0.01 

0.99 
1.25±0.01 
0.85± 0.01 
1.04± 0.05 

1.19 
0.89± 0.01 
0.89± 0.01 

0.88 
0.77± 0.01 
1.33±0.02 
0.73± 0.01 
0.68± 0.05 
0.82± 0.03 
0.52± 0.01 
0.91±0.01 
0.80± 0.01 

0.58 
0.79± 0.01 
0.71±0.01 
0.77± 0.03 
0.61±0.04 
0.62± 0.01 
0.80± 0.02 
0.58± 0.01 
0.59± 0.01 
0.44± 0.01 
0.79± 0.01 
0.75± 0.01 
0.69± 0.01 

0.71 
0.70± 0.01 
0.72± 0.03 
1.28± 0.02 

0.80 
0.71 ±0.02 
0.67± 0.08 
0.55± 0.02 
0.40± 0.05 
1.22± 0.01 

0.92 
0.92± 0.01 

1.09 
0.84± 0.01 

1.11 



Li 

(Lp) 

17.6 
30.4± 0.43 
10.2± 1.65 
23.5± 0.09 
11. 0± 0.18 
18.0±3.00 
8.95± 0.99 

3.36 
13.5± 0.20 

8.30 
5.90± 0.11 
9.58± 0.21 
6.93± 0.04 

8.27 
5.91 ±0.06 

10.27 
6.19± 0.61 

7.08 
3.76± 0.02 

4.59 

3.84 
11.2± 0.10 
8.91±0.62 
3.47± 0.09 
2.78± 0.50 

2.80 
2.70± 0.04 
2.64± 0.04 

2.59 

2.78 
4.65± 0.04 
3.31±0.42 
4.10± 0.60 
2.87± 0.28 
2.84± 0.08 
2.18± 0.08 
2.55± 0.03 
2.23± 0.05 

2.50 

1.37 
3.00± 0.44 
1.78± 0.03 

2.62 

3.12 
1.36± 0.19 
1.82± 0.38 

1.30 
1.97± 0.25 
2.48± 0.41 

3.94 
0.98± 0.16 
5.06± 0.13 

5.65 

7.05 

3.08 
1.1 2± 0.02 

2.39 



L 2 

(Lq) 
4.48± 1.73 
6.14± 0.41 
2.19± 0.35 
5.09± 0.04 
3.71±0.07 
6.00± 1.00 
3.47± 0.38 

0.73 
2.43± 0.07 

0.95 
1.01±0.06 
2.50± 0.09 
2.07± 0.02 

1.76 
3.16± 0.06 

1.34 
2.60± 0.27 

2.49 
0.98± 0.01 
1.01±0.03 

1.28 
0.92± 0.01 
2.88± 0.02 
0.66± 0.03 
0.64± 0.11 
0.91±0.08 
0.42± 0.01 
1.17±0.07 
0.82± 0.03 

0.47 
1.25±0.01 
0.74± 0.02 
0.87± 0.22 
0.48± 0.06 
0.49± 0.02 
0.93± 0.06 
0.41 ±0.03 
0.39± 0.02 
0.25± 0.02 
0.65± 0.02 
0.69± 0.10 
0.47± 0.01 

0.65 
0.74± 0.09 
0.59± 0.05 
1.55±0.32 

0.66 
0.66± 0.04 
0.55±0.10 

0.69 
0.16± 0.03 
2.67± 0.06 

1.59 
1.75±0.07 

1.84 
0.56± 0.01 

1.40 



Mi 
(M e ) 

2.80± 0.60 
2.49± 0.06 
2.09± 0.43 
2.07± 0.05 
1.99± 0.04 
1.94± 0.04 
1.93±0.21 

1.81 
1.80± 0.02 
1.79± 0.14 
1.75±0.04 
1.74± 0.06 
1.74± 0.05 

1.74 
1.73±0.02 
1.72± 0.04 
1.70± 0.16 

1.68 
1.68± 0.03 
1.68± 0.05 

1.64 
1.61±0.05 
1.60± 0.02 
1.57± 0.03 
1.47± 0.15 
1.46± 0.04 
1.45±0.03 
1.43±0.03 
1.43±0.03 

1.42 
1.40± 0.03 
1.38± 0.01 
1.37± 0.06 
1.37± 0.16 
1.35±0.06 
1.31 ±0.05 
1.28± 0.06 
1.27± 0.06 
1.25±0.04 
1.22± 0.03 
1.21±0.15 
1.19± 0.01 

1.18 
1.18± 0.17 
1.06± 0.06 
1.05±0.02 

1.01 
1.00± 0.03 
0.94± 0.14 
0.93± 0.34 
0.72± 0.05 
1.85±0.02 

1.65 
1.64± 0.05 

1.61 
1.59± 0.02 

1.56 



M 2 

(M Q ) 
0.50± 0.17 
0.76± 0.03 
0.45± 0.10 
0.51 ±0.03 
0.60± 0.02 
0.59± 0.02 
0.69± 0.08 

0.29 
0.34± 0.01 
0.14± 0.01 
0.24± 0.02 
0.68± 0.04 
0.59± 0.03 

0.53 
0.85± 0.02 
0.22± 0.02 
0.64± 0.06 

0.49 
0.39± 0.02 
0.28± 0.01 

0.50 
0.15±0.02 
0.51±0.01 
0.31 ±0.02 
0.25± 0.03 
0.43± 0.02 
0.16± 0.01 
0.53± 0.02 
0.35± 0.01 

0.15 
0.33± 0.01 
0.26± 0.01 
0.25± 0.02 
0.17± 0.06 
0.18±0.02 
0.48± 0.03 
0.18± 0.01 
0.17±0.03 
0.10± 0.01 
0.54± 0.02 
0.22± 0.03 
0.29± 0.01 

0.29 
0.24± 0.04 
0.35± 0.03 
0.88± 0.02 

0.52 
0.30± 0.01 
0.13± 0.02 
0.15± 0.05 
0.15±0.04 
0.82± 0.02 

0.30 
0.33± 0.02 

0.82 
0.54± 0.01 

0.85 



M 2 ,i 
(M Q ) 
2.54± 0.69 
2.51±0.13 
2.21±0.38 
2.61± 0.09 
2.35± 0.07 
2.64± 0.15 
2.22± 0.32 

1.86 
2.35± 0.04 
2.1 1± 0.03 
2.05± 0.08 
2.03± 0.15 
2.03±0.10 

2.01 
1.95±0.08 
2.18± 0.06 
2.10± 0.25 

2.24 
1.89± 0.07 
2.01±0.04 

1.88 
2.09± 0.06 
2.29± 0.03 
1.80± 0.08 
1.85±0.16 
1.80± 0.11 
1.79± 0.04 
1.80± 0.10 
1.85±0.05 

1.84 
2.06± 0.03 
1.90± 0.04 
1.98± 0.16 
1.83±0.23 
1.83±0.08 
1.75±0.13 
1.76± 0.06 
1.75±0.11 
1.68± 0.05 
1.47± 0.09 
1.91±0.15 
1.68± 0.04 

1.82 
1.92± 0.17 
1.71±0.13 
1.30± 0.33 

1.51 
1.81 ±0.06 
1.91±0.12 
1.98± 0.15 
1.47± 0.19 
1.88± 0.09 

2.19 
2.21±0.08 

1.60 
1.39± 0.05 

1.29 



M 1A 
(M e ) 

2.11± 0.57 
1.90± 0.10 
1.50± 0.25 
1.37±0.05 
1.40± 0.04 
1.25± 0.07 
1.42± 0.20 

1.28 
1.12± 0.02 
1.13±0.02 
1.1 4± 0.05 
1.29± 0.10 
1.26± 0.06 

1.24 
1.36± 0.06 
1.06± 0.03 
1.21± 0.14 

1.09 
1.1 8± 0.04 
1.1 0± 0.02 

1.18 
0.96± 0.03 
1.00± 0.01 
1.07± 0.05 
0.93± 0.08 
1.00± 0.06 
0.90± 0.02 
1.00± 0.05 
0.93± 0.02 

0.85 
0.82± 0.01 
0.83± 0.02 
0.79± 0.07 
0.81±0.10 
0.80± 0.03 
0.88± 0.07 
0.75± 0.02 
0.74± 0.05 
0.72± 0.02 
0.91± 0.06 
0.64± 0.05 
0.72± 0.02 

0.67 
0.61±0.05 
0.60± 0.05 
0.91± 0.23 

0.68 
0.49± 0.02 
0.34± 0.02 
0.32± 0.02 
0.28± 0.04 
1.50± 0.07 

1.01 
1.01± 0.03 

1.35 
1.30± 0.05 

1.41 



0.83 
0.76 
0.68 
0.52 
0.60 
0.47 
0.64 
0.69 
0.48 
0.53 
0.56 
0.63 
0.62 
0.62 
0.70 
0.49 
0.58 
0.49 
0.62 
0.55 
0.62 
0.46 
0.44 
0.59 
0.50 
0.55 
0.50 
0.56 
0.50 
0.46 
0.40 
0.44 
0.40 
0.45 
0.44 
0.51 
0.43 
0.42 
0.43 
0.62 
0.33 
0.43 
0.37 
0.32 
0.35 
0.70 
0.45 
0.27 
0.18 
0.16 
0.19 
0.80 
0.46 
0.46 
0.84 
0.94 
1.09 



Gy_ 

0.5 

0.5 

0.7 

0.5 

0.6 

0.4 

0.7 

1.2 

0.6 

0.8 

0.9 

0.9 

0.9 

1.0 

1.1 

0.8 

0.9 

0.7 

1.2 

1.0 

1.2 

0.9 

0.7 

1.4 

1.3 

1.4 

1.4 

1.4 

1.3 

1.3 

0.9 

1.2 

1.0 

1.3 

1.3 

1.5 

1.5 

1.5 

1.7 

2.6 

1.1 

1.7 

1.3 

1.1 

1.6 

3.9 

2.4 

1.3 

1.1 

1.0 

2.6 

1.2 

0.8 

0.7 

2.0 

3.1 

4.0 



Ref 

637 

45 

14 

8 

45 
26 

5,13,37 

37 

36 

27,37 
44 
44 
8 

37 
6 

9,37 
40 
37 
8 

35,37 
37 
8 
1 

45 
39 

21,37 

45 

7 

28,37 

37 

24 

12 

25 

29 

7 

45 
44 
8 

37,43 

28,37 

41 

8 

11 

31 

33 

18 

37 

33 

15 

2,37 
14 

3 

37 

37,43 

37 

8 

37 
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Star 


t P 
d 


Ti 
(K) 


Ti 
(K) 


Ri 
(Re) 




R2 

(R©) 






Li 
(L©) 




L 2 
(L©) 




Ml 
(M©) 




M 2 
CM©) 


Af 2 ,i 
(Mq) 






Afi.i 
(M©) 






Ref 


YY Eri 


W 0.321 


5370± 250 


5623± 260 


1.20± 0.06 





77±0 


04 


1 


07± 0.17 





52±0 


08 


1.55±0 


14 


0.62± 0.07 


1.16=1=0 


45 


1 


37± 0.53 


1.18 


5.7 


14 


EROri 


W 0.423 


6200 


6314 


1.39 




1.14 






2.56 




1.85 




1.53 




0.98 


1.05 






1.51 


1.44 


8.1 


10,37 


V502 Oph 


W 0.453 


5888=b 271 


6166± 284 


1.51± 0.14 





93±0 


06 


2 


45± 0.40 


1 


12±0 


IS 


1.51=1= 


28 


0.56== 0.08 


1.73=1=0 


36 


1 


12± 0.23 


0.64 


1.5 


14 


V870 Ara 


W 0.400 


5860± 35 


6210± 35 


1.67± 0.01 





61±0 


01 


2 


96± 0.30 





50=1=0 


01 


1.50=1=0 


01 


0.12± 0.01 


1.89=1=0 


04 





91 ± 0.02 


0.48 


1.2 


33 


BB Peg 


W 0.362 


5780± 7 


6100± 7 


1.28±0.01 





81±0 


01 


1 


61±0.03 





81±0 


01 


1.42=1=0 


02 


0.55± 0.01 


1.57=1=0 


04 


1 


08± 0.03 


0.68 


2.1 


44 


VY Sex 


W 0.443 


5756± 15 


5960± 15 


1.50± 0.01 





86±0 


01 


2 


17± 0.02 





83±0 


01 


1.42=1=0 


02 


0.45± 0.01 


1.73=1=0 


04 





99± 0.02 


0.57 


1.5 


8 


V417 Aql 


W 0.370 


6030 


6256 


1.31 




0.84 






2.02 




0.96 




1.40 




0.50 


1.74 






0.98 


0.57 


1.5 


37 


AEPhe 


W 0.362 


5880± 270 


6160± 280 


1.29± 0.03 





81±0 


02 


1 


74± 0.28 





83±0 


13 


1.38=1=0 


06 


0.63± 0.02 


1.44=1=0 


18 


1 


11± 0.14 


0.77 


2.8 


14 


EZHya 


W 0.450 


5721 


6100 


1.54 




0.85 






2.28 




0.90 




1.37 




0.35 


1.89 






0.85 


0.45 


1.2 


37 


UVLyn 


W 0.415 


6045± 12 


6262± 12 


1.39± 0.01 





89±0 


01 




2.30 


1 


09=1=0 


0? 


1.36=1=0 


02 


0.50== 0.01 


1.80=1=0 


05 





92± 0.02 


0.51 


1.4 


35,37 


AH Vir 


W 0.407 


5300 


5671 


1.41 




0.84 






1.40 




0.65 




1.36 




0.41 


1.67 






0.94 


0.56 


1.7 


37 


V842 Her 


W 0.419 


6000 


6280 


1.46 




0.81 






2.47 




0.92 




1.36 




0.35 


1.90 






0.84 


0.44 


1.2 


37 


WUMa 


W 0.334 


6026± 278 


6310± 291 


1.17± 0.03 





85±0 


02 


1 


62± 0.26 


1 


02±0 


17 


1.35=1=0 


09 


0.69± 0.05 


1.45=1=0 


32 


1 


09± 0.24 


0.75 


2.7 


14 


QWGem 


W 0.358 


5890± 10 


6100± 10 


1.26± 0.03 





75±0 


02 




1.71 





69=1=0 


04 


1.31=1=0 


04 


0.44± 0.01 


1.66=1=0 


06 





90± 0.03 


0.54 


1.8 


21,37 


SS Ari 


W 0.406 


5860 


6123 


1.36 




0.80 






1.96 




0.81 




1.31 




0.40 


1.78 






0.84 


0.47 


1.4 


37 


V752 Cen 


W 0.370 


5955 


6221 


1.28 




0.75 






1.85 




0.76 




1.30 




0.40 


1.76 






0.84 


0.48 


1.5 


4,37 


AM Leo 


W 0.366 


6100± 5 


6221± 5 


1.23± 0.01 





85±0 


01 


1 


84± 0.02 





95±0 


01 


1.29=1=0 


01 


0.59± 0.01 


1.59=1=0 


04 





95 ± 0.03 


0.60 


2.0 


45 


V781 Tau 


W 0.408 


5885± 50 


6150± 50 


1.21± 0.01 





85±0 


01 


1 


58± 0.04 





93±0 


04 


1.29=1=0 


07 


0.57± 0.03 


1.62=1=0 


13 





94± 0.08 


0.58 


1.9 


18,38 


VI 191 Cyg 


W0.313 


6500 


6610 


1.31 




0.52 






2.71 




0.46 




1.29=1=0 


08 


0.13± 0.01 


1.85=1=0 


03 





71±0.01 


0.38 


1.3 


34 


RZCom 


W 0.339 


6166± 284 


6457± 298 


1.12± 0.03 





78±0 


02 


1 


62± 0.26 





93±0 


15 


1.23=1=0 


09 


0.55± 0.04 


1.65=1=0 


23 





86± 0.12 


0.52 


1.8 


14 


GM Dra 


W 0.339 


6306± 58 


6450± 58 


1.25± 0.01 





61±0 


01 


2 


19± 0.17 





56±0 


01 


1.21=1=0 


04 


0.22± 0.02 


1.83=1=0 


07 





67± 0.02 


0.36 


1.3 


7 


FUDra 


W 0.307 


5800± 8 


6133± 8 


1.12± 0.01 





61±0 


01 




1.27 





46=1=0 


02 


1.17=1=0 


04 


0.29± 0.01 


1.68=1=0 


05 





70± 0.02 


0.42 


1.7 


35,37 


UPeg 


W 0.375 


5860± 10 


5841± 10 


1.22± 0.01 





74±0 


01 




1.58 





58=1=0 


02 


1.15=1=0 


01 


0.38±0.01 


1.66=1=0 


05 





72± 0.02 


0.43 


1.8 


28,37 


GZAnd 


W 0.305 


5810± 25 


6200± 25 


LOliO.Ol 





74±0 


01 




1.03 





73±0 


03 


1.12=1=0 


02 


0.59± 0.02 


1.44=1=0 


10 





83± 0.06 


0.58 


2.8 


3,37 


AO Cam 


W 0.330 


5590± 33 


5900± 33 


1.09± 0.01 





73±0 


01 


1 


03± 0.07 





58±0 


01 


1.12=1=0 


01 


0.49± 0.01 


1.50=1=0 


05 





78± 0.02 


0.52 


2.5 


3.7 


TWCet 


W 0.317 


5450 


5630 


0.99 




0.76 






0.77 




0.52 




1.06 




0.61 


1.18 






0.87 


0.73 


5.3 


37 


BVDra 


W 0.350 


6245± 30 


6345 ± 30 


1.12± 0.01 





76±0 


01 




1.71 





84±0 


04 


1.04=1= 


02 


0.43± 0.01 


1.76=1=0 


05 





59± 0.02 


0.34 


1.5 


19,37 


BX Peg 


W 0.280 


5300± 10 


5528± 10 


0.97± 0.01 





62±0 


01 





66± 0.01 





32±0 


01 


1.02=1=0 


12 


0.38± 0.05 


1.40=1=0 


20 





68±0.10 


0.48 


3.1 


32 


OU Ser 


W 0.297 


5960=1= 42 


6380± 42 


1.09± 0.01 





51±0 


01 




1.34 





38±0 


02 


1.02=1=0 


01 


0.18±0.01 


1.73=1=0 


05 





50± 0.02 


0.29 


1.6 


27,37 


AB And 


W 0.332 


5495± 250 


5888± 270 


1.05± 0.01 





76±0 


01 





87± 0.20 





59±0 


12 


1.01=±= 


02 


0.49± 0.01 


1.50=1=0 


14 





67± 0.06 


0.45 


2.4 


16 


SWLac 


W 0.321 


5347 


5630 


1.03 




0.94 






0.78 




0.80 




0.98 




0.78 


0.94 






0.93 


0.99 


11.7 


37 


TYBoo 


W 0.317 


5800± 200 


6180±200 


1.00± 0.01 





69±0 


01 


1 


02± 0.14 





62±0 


09 


0.93=1= 


02 


0.40=1= 0.01 


1.66=1=0 


10 





51± 0.03 


0.30 


1.8 


30 


VW Cep 


W 0.278 


5012± 231 


5248± 242 


0.91±0.02 





62±0 


01 





47± 0.08 





26±0 


04 


0.93=1= 


02 


0.40± 0.01 


1.27=1=0 


11 





64± 0.05 


0.50 


4.2 


14 


BW Dra 


W 0.292 


5980=1= 10 


6164± 10 


0.98± 0.01 





55±0 


01 




1.10 





39±0 


02 


0.92=1= 


02 


0.26± 0.01 


1.65=1=0 


05 





45± 0.01 


0.27 


1.8 


19,37 


TXCnc 


W 0.383 


5888=b 271 


6165± 284 


1.12± 0.03 





83±0 


04 


1 


35± 0.22 





91±0 


15 


0.91=1=0 


10 


0.50± 0.06 


1.71=1=0 


29 





50± 0.08 


0.29 


1.6 


14 


V757 Cen 


W 0.343 


5900 


6000 


0.97 




0.80 






1.02 




0.74 




0.88 




0.59 


1.45 






0.59 


0.41 


2.7 


37 


V829 Her 


W 0.358 


5900± 21 


5380± 21 


1.07± 0.01 





74±0 


01 




1.25 





42±0 


02 


0.86=1= 


02 


0.37=1= 0.01 


1.53=1=0 


05 





47± 0.02 


0.31 


2.3 


37,43 


CC Com 


W 0.221 


4200± 180 


4300± 180 


0.71±0.01 





53±0 


01 





14± 0.01 





09±0 


01 


0.72=1= 


02 


0.38=1= 0.01 


0.86=1= 


11 





56± 0.07 


0.65 


15.9 


20 


XY Leo 


W 0.284 


4361± 8 


4800± 8 


0.83± 0.13 





66±0 


10 





22± 0.08 





21±0 


07 


0.76=1= 


15 


0.46± 0.06 


1.03=1=0 


46 





57± 0.26 


0.55 


8.5 


33 


V523 Cas 


W 0.234 


4410± 5 


4736± 5 


0.74± 0.04 





55±0 


02 





18± 0.02 





13±0 


01 


0.75=1= 


03 


0.38± 0.02 


1.00±0 


13 





54± 0.07 


0.55 


9.6 


42 


BH Cas 


W 0.406 


4790=b 100 


4980± 100 


1.05± 0.01 





75±0 


01 





52± 0.01 





31±0 


01 


0.74=1= 


06 


0.35± 0.03 


1.43=1=0 


12 





38± 0.03 


0.26 


2.8 


22,37 


RW Dor 


W 0.285 


4780 


5200 


0.80 




0.67 






0.30 




0.30 




0.64 




0.43 


1.28 






0.35 


0.28 


4.1 


37 


RWCom 


W 0.237 


5120± 10 


5400± 10 


0.71± 0.01 





46±0 


01 




0.31 





16±0 


01 


0.56=1= 


06 


0.20=1= 0.03 


1.41=1=0 


12 





15± 0.01 


0.11 


3.0 


23,37 
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